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Abstract 

This is a short review in honor of B. Mandelbrot's 80st birthday, to 
appear in Wilmott magazine. We discuss how multiplicative cascades 
and related multifractal ideas might be relevant to model the main 
statistical features of financial time series, in particular the intermit- 
tent, long-memory nature of the volatility. We describe in details the 
Bacry-Muzy-Delour multifractal random walk. We point out some 
inadequacies of the current models, in particular concerning time re- 
versal symmetry, and propose an alternative family of multi-timescale 
models, intermediate between GARCH models and multifractal mod- 
els, that seem quite promising. 
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1 Introduction 



Financial time series represent an extremely rich and fascinating source of ques- 
tions, for they store a quantitative trace of human activity, sometimes over hun- 
dreds of years. These time series, perhaps surprisingly, turn out to reveal a very 
rich and particularly non trivial statistical texture, like fat tails and intermittent 
volatility bursts, have been by now well-characterized by an uncountable number 
of empirical studies, initiated by Benoit Mandelbrot more than fourty years ago 
PP. Quite interestingly, these statistical anomalies are to some degree universal, 
i.e, common across different assets (stocks, currencies, commodities, etc.), regions 
(U.S., Europe, Asia) and epochs (from wheat prices in the 18th century to oil 
prices in 2004). The statistics of price changes is very far from the Bachelier- 
Black-Scholes random walk model which, nevertheless, is still today the central 
pillar of mathematical finance: the vast majority of books on option pricing writ- 
ten in the last ten years mostly focus on Brownian motion models and up until 
recently very few venture into the wonderful world of non Gaussian random walks 
[2j. This is the world, full of bushes, gems and monsters, that Mandelbrot started 
exploring for us in the sixties, charting out its scrubby paths, on which droves of 
scientists - admittedly with some delay - now happily look for fractal butterflies, 
heavy-tailed marsupials or long- memory elephants. One can only be superlative 
about his relentless efforts to look at the world with new goggles, and to offer the 
tools that he forged to so many different communities: mathematicians, physicists, 
geologists, meteorologists, fractologists, economists - and, for our purpose, finan- 
cial engineers. In our view, one of Mandelbrot's most important methodological 
messages is that one should look at data, charts and graphs in order to build one's 
intuition, rather than trust blindly the result of statistical tests, often inadequate 
and misleading, in particular in the presence of non Gaussian effects [3]. This 
visual protocol is particularly relevant when modeling financial time series: as dis- 
cussed below, well chosen graphs often allow one to identify important effects, rule 
out an idea or construct a model. This might appear as a trivial statement but, 
as testified by Mandelbrot himself, is not: he fought all his life against the Bour- 
baki principle that pictures and graphs betray [3]. Unreadable tables of numbers, 
flooding econometrics papers, perfectly illustrate that figures remain suspicious in 
many quarters. 

The aim of this paper is to give a short review of the stylized facts of financial 
time series, and of how multifractal stochastic volatility models, inherited from 
Kolmogorov's and Mandelbrot's work in the context of turbulence, fare quite well 
at reproducing many important statistical features of price changes. We then 
discuss some inadequacies of the current models, in particular concerning time 
reversal symmetry, and propose an alternative family of multi-timescale models, 
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intermediate between GARCH models and multifractal models, that seem quite 
promising. We end on a few words about how these models can be used in practice 
for volatility filtering and option pricing. 



2 Universal features of price changes 



The modeling of random fluctuations of asset prices is of obvious importance in 
finance, with many applications to risk control, derivative pricing and systematic 
trading. During the last decade, the availability of high frequency time series has 
promoted intensive statistical studies that lead to invalidate the classic and popular 
"Brownian walk" model, as anticipated by Mandelbrot [3]. In this section, we 
briefly review the main statistical properties of asset prices that can be considered 
as universal, in the sense that they are common across many markets and epochs 



If one denotes p(t) the price of an asset at time t, the return r T (t), at time t 
and scale r is simply the relative variation of the price from t to t + r: r T (t) = 
\p(t + r) - p(t)} /p{t) ~ hxp(t + r) - lnp(t). 

The simplest "universal" feature of financial time series, uncovered by Bachelier 
in 1900, is the linear growth of the variance of the return fluctuations with time 
scale. More precisely, if mr is the mean return at scale r, the following property 
holds to a good approximation: 



where (...) e denotes the empirical average. This behaviour typically holds for r 
between a few minutes to a few years, and is equivalent to the statement that 
relative price changes are, to a good approximation, uncorrelated. Very long time 
scales (beyond a few years) are difficult to investigate, in particular because the 
average drift m becomes itself time dependent, but there are systematic effects 
suggesting some degree of mean-reversion on these long time scales. 

The volatility a is the simplest quantity that measures the amplitude of return 
fluctuations and therefore that quantifies the risk associated with some given as- 
set. Linear growth of the variance of the fluctuations with time is typical of the 
Brownian motion, and, as mentioned above, was proposed as a model of market 
fluctuations by Bachelier. 1 In this model, returns are not only uncorrelated, as 
mentioned above, but actually independent and identical Gaussian random vari- 
ables. However, this model completely fails to capture other statistical features 

In Bachelier's model, absolute price changes, rather than relative returns, were con- 
sidered; there are however only minor differences between the two at short time scales, 
t < 1 month, but see the detailed discussion [7], Ch. 7. 
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Figure 1 : 1-a: Absolute value of the daily price returns for the Dow- Jones index 
over a century (1900-2000), and zoom on different scales (1990-2000 and 1995). 
Note that the volatility can remain high for a few years (like in the early 1930's) 
or for a few days. This volatility clustering can also be observed on high frequency 
(intra-day) data. 1-b: Same plot for a Brownian random walk, which shows a 
featureless pattern in this case. 

of financial markets that even a rough analysis of empirical data allows one to 
identify, at least qualitatively: 

(i) The distribution of returns is in fact strongly non Gaussian and its shape 
continuously depends on the return period r: for r large enough (around few 
months), one observes quasi-Gaussian distributions while for small r values, 
the return distributions have a strong kurtosis (see Fig. EJ). Several care- 
ful studies suggest that these distributions can be characterized by Pareto 
(power-law) tails Irl" 1 "^ with an exponent /i in the range 3 — 5 even for liq- 
uid markets such as the US stock index, major currencies, or interest rates 
[HI EH El DHj • Note that fj, > 2, and excludes an early suggestion by Mandel- 
brot that security prices could be described by Levy stable laws. Emerging 
markets, however, have even more extreme tails, with an exponent \i that 
can be less than 2 - in which case the volatility is formally infinite - as found 
by Mandelbrot in his famous study of cotton prices A natural conjec- 
ture is that as markets become more liquid, the value of fj, tends to increase 
(although illiquid, yet very regulated markets, could on the contrary show 
truncated tails because of artificial trading rules |llj). 

(ii) Another striking feature is the intermittent and correlated nature of return 
amplitudes. At some given time period r, the volatility is a quantity that 
can be defined locally in various ways: the simplest ones being the square 
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return r^.(t) or the absolute return |r r (i)|, or a local moving average of these 
quantities. The volatility signals are characterized by self-similar outbursts 
(see Fig. ^) that are similar to intermittent variations of dissipation rate in 
fully developed turbulence • The occurrence of such bursts are strongly 
correlated and high volatility periods tend to persist in time. This feature 
is known as volatility clustering (131 1141 |S1 E] • This effect can analyzed more 
quantitatively: the temporal correlation function of the (e.g. daily) volatility 
can be fitted by an inverse power of the lag, with a rather small exponent 
in the range 0.1-0.3 [1111311111711171. 

(iii) One observes a non-trivial "multifractal" scaling [T%1 H51 12U1 12T1 IT71 1221 1531 • 

in the sense that higher moments of price changes scale anomalously with 
time: 

M q (r) = (\r T (t) - mr\ q ) e ~ A g r<\ (2) 

where the index £ g is not equal to the Brownian walk value q/2. As will be 
discussed more precisely below, this behaviour is intimately related to the 
intermittent nature of the volatility process. 

(iv) Past price changes and future volatilities are negatively correlated, at least 
on stock markets. This is the so called leverage effect, which reflects the fact 
that markets become more active after a price drop, and tend to calm down 
when the price rises. This correlation is most visible on stock indices |24j . 
This leverage effect leads to an anomalous negative skew in the distribution 
of price changes as a function of time. 

The most important message of these empirical studies is that price changes 
behave very differently from the simple geometric Brownian motion: extreme 
events are much more probable, and interesting non linear correlations (volatility- 
volatility and price-volatility) are observed. These "statistical anomalies" are very 
important for a reliable estimation of financial risk and for quantitative option 
pricing and hedging, for which one often requires an accurate model that cap- 
tures return statistical features on different time horizons r. It is rather amazing 
to remark that empirical properties (i-iv) are, to some extent, also observed on 
experimental velocity data in fully developed turbulent flows (see Fig. |2J) - The 
framework of scaling theory and multifractal analysis, initially proposed to char- 
acterize turbulent signals by Mandelbrot and others ^2], may therefore be well- 
suited to further characterize statistical properties of price changes on different 
time periods. 
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Figure 2: Continuous deformation of turbulent velocity increments and financial 
returns distributions from small (top) to large (bottom) scales. 2-a: Standard- 
ized probability distribution functions of spatial velocity increments at different 
length scales in a high Reynolds number wind tunnel turbulence experiment. The 
distributions are plotted in logarithmic scale so that a parabola corresponds to a 
Gaussian distribution. 2-b: Standardized p.d.f. of S&P 500 index returns at dif- 
ferent time scales from few minutes to one month. Notice that because of sample 
size limitation, the noise amplitude is larger than in turbulence. 



3 From multifractal scaling to cascade pro- 
cesses 

3.1 Multi-scaling of asset returns 

For the geometric Brownian motion, or for the Levy stable processes first sug- 
gested by Mandelbrot as an alternative, the return distribution is identical (up to 
a rescaling) for all r. As emphasized in previous section, the distribution of real 
returns is not scale invariant, but rather exhibits multi-scaling. Fig. |3] illustrates 
the empirical multifractal analysis of S&P 500 index return. As one can see in Fig. 
Efb), the scaling behavior © that corresponds to a linear dependence in a log-log 
representation of the absolute moments versus the time scale r, is well verified 
over some range of time scales (typically 3 decades) . 

The multifractal nature of the index fluctuations can directly be checked in 
Fig. Et c )> where one sees that the moment ratios strongly depend on the scale 
r. The estimated Cq spectrum (Fig. Efd)) has a concave shape that is well fitted 
by the parabola: ( g = q(l/2 + A 2 ) - X 2 q 2 /2 with A 2 ~ 0.03. The coefficient 



6 



A 2 that quantifies the curvature of the C g function is called, in the framework of 
turbulence theory, the intermittency coefficient. The most natural way to account 
for the multi-scaling property @ is through the notion of cascade from coarse to 
fine scales. 

3.2 The cascade picture 

As noted previously, for the geometric Brownian motion, the return probability 
distributions at different time periods r are Gaussian and thus differ only by their 
width that is proportional to \fr. If x(t) is a Brownian motion, this property can 
be shortly written as 

r T (t) =iaw o- T e(t) (3) 
o- ST = s 1/2 a T (4) 

where =i aw means that the two quantities have the same probability distribution, 
e(t) is a standardized Gaussian white noise and a T is the volatility at scale r. When 
going from some scale r to the scale st, the return volatility is simply multiplied by 
s 1 / 2 . The cascade model assumes such a multiplicative rule but the multiplicative 
factor is now a random variable and the volatility itself becomes a random process 
a T (t): 

r T (t) =i aw cr T (t)e(t) (5) 

(Tsrit) =law W s O T (t) (6) 

where the law of W s depends only on the scale ratio s and is independent of o~ T (t). 
Let T be some coarse time period and let s < 1. Then, by setting W = et, and by 
iterating equation (JBJ) n times, one obtains: 

a s n T {t) = W s n(j T {t) = e^ n ^ il ^a T {t) (7) 

law law 

Therefore, the logarithm of the volatility at some fixed scale r = s n T, can be 
written as a sum of an arbitrarily large number n of independent, identically 
distributed random variables. Mathematically, this means that the logarithm of 
the volatility (and hence £, the logarithm of the "weights" W) belongs the the 
class of the so-called infinitely divisible distributions [2^1 . The simplest of such 
distributions (often invoked using the central limit theorem) is the Gaussian law. 
In that case, the volatility is a log-normal random variable. As explained in the 
next subsection, this is precisely the model introduced by Kolmogorov in 1962 
to account for the intermittency of turbulence |26| . It can be proven that the 
random cascade equations © and © directly lead to the deformation of return 
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probability distribution functions as observed in Fig. Using the fact that £ s is 
a Gaussian variable of mean /zln(s) and variance A 2 ln(s), one can compute the 
absolute moment of order q of the returns at scale r = s n T (s < 1). One finds: 

(\r T (t)\«)=A q (-) (8) 

where A q = (\r T (t)\ q ). 

Using a simple multiplicative cascade, we thus have recovered the empirical 
findings of previous sections. For log-normal random weights W s , the return pro- 
cess is multifractal with a parabolic £ g scaling spectrum: Q = q\i — X 2 q 2 /2 where 
the parameter \i is related to the mean of In W s and the curvature A 2 (the intermit- 
tency coefficient) is related to the variance of lnPU s and therefore to the variance 
of the log- volatility ln(<7 T ): 

((lna r (t)) 2 ) - (lna T (t)} 2 = -A 2 ln(r) + V (9) 

The random character of In W s is therefore directly related to the intermittency of 
the returns. 

4 The multifractal random walk 

The above cascade picture assumes that the volatility can be constructed as a 
product of random variables associated with different time scales. In the previous 
section, we exclusively focused on return probability distribution functions (mono- 
variate laws) and scaling laws associated with such models. Explicit constructions 
of stochastic processes which marginals satisfy a random multiplicative rule, were 
first introduced by Mandelbrot .27, and are known as Mandelbrot's cascades or 
random multiplicative cascades. The construction of a Mandelbrot cascade, illus- 
trated in Fig. 01 always involves a discrete scale ratio s (generally s = 1/2). One 
begins with the volatility at the coarsest scale and proceeds in a recursive manner 
to finer resolutions: The n-th step of the volatility construction corresponds to 
scale 2~ n and is obtained from the (n — l)-th step by multiplication with a pos- 
itive random process W the law of which does not depend on n. More precisely, 
the volatility in each subinterval is the volatility of its parent interval multiplied 
by an independent copy of W. 

Mandelbrot cascades are considered to be the paradigm of multifractal pro- 
cesses and have been extensively used for modeling scale-invariance properties in 
many fields, in particular statistical finance by Mandelbrot himself I22j. How- 
ever, this class of models presents several drawbacks: (i) they involve a preferred 
scale ratio s, (ii) they are not stationary and (iii) they violate causality. In that 
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respect, it is difficult to see how such models could arise from a realistic (agent 
based) description of financial markets. 

Recently, Bacry, Muzy and Delour (BMD) |17l I29j introduced a model that 
does not possess any of the above limitations and captures the essence of cascades 
through their correlation structure (see also |311 I32| for alternative multifractal 
models). 

In the BMD model, the key ingredient is the volatility correlation shape that 
mimics cascade features. Indeed, as remarked in ref. [US], the tree like structure 
underlying a Mandelbrot cascade, implies that the volatility logarithm covariance 
decreases very slowly, as a logarithm function, i.e., 

(Ha T (t)) ln(<r T (t + At))) - (ln(a T (t))) 2 = C - A 2 ln(At + r) (10) 

This equation can be seen as a generalization of Kolmogorov equation that 
describes only the behavior of the log-volatility variance. It is important to note 
that such a logarithmic behaviour of the covariance has indeed been observed for 
empirically estimated log- volatilities in various stock market data [SSI- The BMD 
model involves Eq. ()10|) within the continuous time limit of a discrete stochastic 
volatility model. One first discretizes time in units of an elementary time step tq 
and sets t = itq. The volatility Oi at "time" i is a log- normal random variable 
such that <7j = (Toexp^j, where the Gaussian process & has the same covariance 
as in Eq. (fTU|) : 

(&) = -A 2 ln(-) = Mo ; m) - ill = A 2 ln(-) - A 2 - j\ + 1), (11) 

TO T 

for \i — ]\tq < T. Here T is a large cut-off time scale beyond which the volatility 
correlation vanishes. In the above equation, the brackets stands for the mathe- 
matical expectation. The choice of the mean value /io is such that (a 2 ) = <Tq. As 
before, the parameter A 2 measures the intensity of volatility fluctuations (or, in 
the finance parlance, the 'vol of the vol'), and corresponds to the intermittency 
parameter. 

Now, the price returns are constructed as: 

x [(i + l)r ] - x [itq] = r T0 (i) = OiEi = a e^£i, (12) 

where the are a set of independent, identically distributed random variables of 
zero mean and variance equal to To- One also assumes that the and the £j are 
independent (but see [31] where some correlations are introduced). In the original 
BMD model, Ej's are Gaussian, and the continuous time limit tq = dt — > is 
taken. Since x = hip, where p is the price, the exponential of a sample path of 
the BMD model is plotted in Fig. |SJ(a), which can be compared to the real price 
chart ofEta). 
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The multifractal scaling properties of this model can be computed explicitly. 
Moreover, using the properties of multivariate Gaussian variables, one can get 
closed expressions for all even moments M q {r) {q = 2k). In the case q = 2 one 
trivially finds: 

M 2 (r = £t ) = a 2 £r = a 2 r, (13) 

independently of A 2 . For q ^ 2, one has to distinguish between the cases qX 2 < 1 
and q\ 2 > 1. For qX 2 < 1, the corresponding moments are finite, and one finds, 
in the scaling region To <C r < T, a genuine multifractal behaviour |171 I29| . for 
which Eq. © holds exactly. For qX 2 > 1, on the other hand, the moments diverge, 
suggesting that the unconditional distribution of x(t + r) — x(t) has power-law tails 
with an exponent \i = 1/X 2 (possibly multiplied by some slow function). These 
multifractal scaling properties of BMD processes are numerically checked in Figs. 
|5J-a and[5J-b where one recovers the same features as in Fig. |3]for the S&P 500 index. 
Since volatility correlations are absent for t > T, the scaling becomes that of a 
standard random walk, for which £ 9 = q/2. The corresponding distribution of price 
returns thus becomes progressively Gaussian. An illustration of the progressive 
deformation of the distributions as r increases, in the BMD model is reported 
in Fig. [SJ-c. This figure can be directly compared to Fig. [2J As shown in Fig. 
03-e, this model also reproduces the empirically observed logarithmic covariance of 
log- volatility (Eq. (|1U|) . This functional form, introduced at a "microscopic" level 
(at scale To — ► 0), is therefore stable across finite time scales r. 

To summarize, the BMD process is attractive for modeling financial time series 
since it reproduces most "stylized facts" reviewed in section 2 and has exact mul- 
tifractal properties as described in section 3. Moreover, this model has stationary 
increments, does not exhibit any particular scale ratio and can be formulated in 
a purely causal way: the log volatility £j can be expressed as a sum over "past" 
random shocks, with a memory kernel that decays as the inverse square root of the 
time lag [351 • Let us mention that generalized stationary continuous cascades with 
compound Poisson and infinitely divisible statistics have been recently contructed 
and mathematically studied by Mandelbrot and Barral |36j . and in P0 However, 
there is no distinction between past and future in this model - but see below. 

5 Multifractal models: empirical data under 
closer scrutiny 

Multifractal models either postulate or predict a number of precise statistical fea- 
tures that can be compared with empirical price series. We review here several 
successes, but also some failures of this family of models, that suggest that the 
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century old quest for a faithful model of financial prices might not be over yet, 
despite Mandelbrot's remarkable efforts and insights. Our last section will be de- 
voted to an idea that may possibly bring us a little closer to the yet unknown 
"final" model. 

5.1 Volatility distribution 

Mandelbrot's cascades construct the local volatility as a product of random volatil- 
ities over different scales. As mentioned above, this therefore leads to the local 
volatility being log-normally distributed (or log infinitely divisible) , an assumption 
that the BMD model postulates from scratch. Several direct studies of the distri- 
bution of the volatility are indeed compatible with log-normality; option traders 
actually often think of volatility changes in relative terms, suggesting that the log- 
volatility is the natural variable. However, other distributions, such as an inverse 
gamma distribution fits the data equally well, or even better |37l 17] . 

5.2 Intermittency coefficient 

One of the predictions of the BMD model is the equality between the intermittency 
coefficient estimated from the curvature of the £ g function and the slope of the log- 
volatility covariance logarithmic decrease. The direct study of various empirical 
log-volatility correlation functions, show that can they can indeed be fitted by 
a logarithm of the time lag, with a slope that is in the same ball-park as the 
corresponding intermittency coefficient A 2 . The agreement is not perfect, though. 
These empirical studies furthermore suggest that the integral time T is on the 
scale of a few years [T7| . 

5.3 Tail index 

As mentioned above, and as realized by Mandelbrot long ago through his famous 
"star equation" [27], multifractal models lead to power law tails for the distribution 
of returns. The BMD model predicts that the power-law index should be [i = 1/A 2 , 
which, for the empirical values of A 2 ~ 0.03, leads to a value of \i ten times larger 
than the empirical value of /i, found in the range 3 — 5 for many assets (HI El lEj • 
Of course, the random variable e could itself be non Gaussian and further fattens 
the tails. However, as recently realized by Bacry, Kozhemyak and Muzy [HE], a 
kind of 'ergodicity breaking' seems to take place in these models, in such a way 
that the theoretical value /i = 1/A 2 may not correspond to the most probable 
value of the estimator of the tail of the return distribution for a given realization 
of the volatility, the latter being much smaller than the former. So, even if this 



11 



scenario is quite non trivial, multifractal models with Gaussian residuals may still 
be consistent with fat tails. 

5.4 Time reversal invar iance 

Both Mandelbrot's cascade models and the BMD model are invariant under time 
reversal symmetry, meaning that no statistical test of any nature can distinguish 
between a time series generated with such models and its time reversed. There 
are however various direct indications that real price changes are not time reversal 
invariant. One well known fact is the "leverage effect" (point (iv) of Section 2, see 
also |24j). whereby a negative price change is on average followed by a volatility 
increase. This effect leads to a skewed distribution of returns, which in turn 
can be read as a skew in option volatility smiles. In the BMD model, one has 
by symmetry that all odd moments of the process vanish. A simple possibility, 
recently investigated in |34| . is to correlate negatively the variable with 'past' 
values of the variables Sj, j < i, through a kernel that decays as a power-law. 
In this case, the multifractality properties of the model are preserved, but the 
expression for £ q is different for q even and for q odd (see also |39j ) . 

Another effect breaking the time reversal invariance is the asymmetric struc- 
ture of the correlations between the past and future volatilities at different time 
horizons. This asymmetry is present in all time series, even when the above lever- 
age effect is very small or inexistent, such as for the exchange rate between two 
similar currencies (e.g. Dollar vs. Euro). The full structure of the volatility corre- 
lations is shown by the following "mug shots" introduced in [301 • These plots show 
how past volatilities measured on different time horizons (horizontal axis) affect 
future volatilities, again on different time horizons (vertical axis). For empirical 
price changes, the correlations between past volatilities and future volatilities is 
asymmetric, as shown on Fig. |HJ This figure reveals clearly two very important 
features: (a) the volatility at a given time horizon is affected mainly by the volatili- 
ties at longer time horizon (the "cascade" effect, partly discovered in [33]), and (b) 
correlations are strongest for time horizons corresponding to the natural human 
cycles, i.e: intra-day, one day, one week and one month. The feature (a) points 
to a difference with the volatility cascade used in turbulence, where eddies at a 
given scale break down into eddies at the immediately smaller scale. For financial 
time series, the volatilities at all time horizons feed the volatility at the shortest 
time horizon [I^l • The overall asymmetry around the diagonal of the figure mea- 
sures the asymmetry of the price time series under time reversal invariance. For a 
time reversal symmetric process, these mug shots are symmetric. An example of 
a very symmetric process is given in Fig. [JJ for a volatility cascade as given in Eq. 
(7). The process for the "log- volatility" £ n at the n-th scale obeys an Ornstein- 
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Uhlenbeck process, with a characteristic time s~ n . In order to obtain a realistic 
cascade of volatilities, the (instantaneous) mean for the process £ n is taken as the 
volatility at the larger scale Yet, the resulting mug shot is very different from 

the empirical one, in particular it is exactly symmetric (at least theoretically). 

Building a multifractal model that reproduces the observed time asymmetry is 
still an open problem. There might however be other avenues, as we now discuss. 

6 Multi-timescale GARCH processes and sta- 
tistical feedback 

The above description of financial data using multifractal, cascade-like ideas is 
still mostly phenomenological. An important theoretical issue is to understand the 
underlying mechanisms that may give rise to the remarkable statistical structure 
of the volatility emphasized by multifractal models: (nearly) log-normal, with 
logarithmically decaying correlations. Furthermore, from a fundamental point 
of view, the existence of two independent statistical processes, one for returns 
and another for the volatility, may not be very natural. In stochastic volatility 
models, such as the multifractal model, the volatility has its own dynamics and 
sources of randomness, without any feedback from the price behaviour. The time 
asymmetry revealed by the above mug-shots however unambiguously shows that 
a strong feedback, beyond the leverage effect, exists in financial markets. 

This time asymmetry, on the other hand, is naturally present in GARCH 
models, where the volatility is a (deterministic) function of the past price changes. 
The underlying intuition is that when recent price changes are large, the activity of 
market agents increases, in turn creating possibly large price changes. This creates 
a non linear feedback, where rare events trigger more rare events, generating non 
trivial probability distributions and correlations. Most GARCH models describe 
the feedback as quadratic in past returns, with the following general form. The 
(normalized) return at time i, computed over a time interval fcro is given by = 
(xi — Xi-k)/y/kro- The general form of the volatility is: 

j<i k>0 

where o~q is a 'bare' volatility that would prevail in the absence of feedback, and 
L(i — j; k) a kernel that describes how square returns on different time scales k, 
computed for different days j in the past, affect the uncertainty of the market at 
time i. [Usually, these processes are rather formulated through a set of recursive 
equations, with a few parameters. By expanding these recursion equations, one 
obtains the AR(oo) form given above.] 
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The literature on GARCH-like processes is huge, with literaly dozen of varia- 
tions that have been proposed (for a short perspective, see @U1I^)- Some of these 
models are very successful at reproducing most of the empirical facts, including the 
"mug shot" for the volatility correlation. An example, that was first investigated 
in |4L)| and recently rediscovered in |42j as a generalization of previous work |43j . 
is to choose L(i — j;k) = 5ijK(k). This choice means that the current volatility 
is only affected by observed price variations between different dates i — k in the 
past and today i. The quadratic dependence of the volatility on past returns is 
precisely the one found in the framework of the statistical feedback process in- 
troduced in [131) w here the volatility is proportional to a negative power of the 
probability of past price changes - hence the increased volatility after rare events. 
The time series resulting from the above choice of kernel are found to exhibit fat, 
power-law tails, volatility bursts and long range memory |42| . Furthermore, the 
distribution of volatilities is found to be very close to log-normal [12] . Choosing K 
to be a power-law of time |4U1 144] allows one to reproduce the power-law decay of 
volatility correlations observed in real data ^2J EE]- the power-law response 
of the volatility to external shocks |351 145j , an apparent multifractal scaling |46| 
and, most importantly for our purpose, the shape of the mug-shots shown above 
|4()j . Many of these results can be obtained analytically |421 144] . Note that one 
could choose K{k) to spike for the day, week, month time scales unveiled by the 
mug-shots [3011^. Yet, a detailed systematic comparison of empirical facts against 
the predictions of the models, GARCH-like or multifractal-like, is still lacking. 

The above model can be generalized further by postulating a Landau-like ex- 
pansion of the volatility as a function of past returns as: 

fc>0 

with G(r) = g\r + giT 1 + ... [II]. The case g\ = corresponds to the model 
discussed above, whereas g\ < allows one to reproduce the leverage effect and 
a negative skewness. Other effects, such as trends, could also be included |17] . 
The interest of this type of models, compared to multifractal stochastic volatility 
models, is that their fundamental justification, in terms of agent based strategy, is 
relatively direct and plausible. One can indeed argue that agents use thresholds for 
their interventions (typically stop losses, entry points or exit points) that depend 
on the actual path of the price over some investment horizon, which may differ 
widely between different investors - hence a power-law like behaviour of the kernels 
K(k). Note in passing that, quite interestingly, Eq. ((11)) is rather at odds with the 
efficient market hypothesis, which asserts that the price past history should have 
no bearing whatsoever on the behaviour of investors. Establishing the validity of 
Eq. (|14jl could have important repercussions on that front, too. 



14 



7 Conclusion: why are these models useful 
anyway ? 

For many applications, such as risk control and option pricing, we need a model 
that allows one to predict, as well as possible, the future volatility of an asset over 
a certain time horizon, or even better, to predict that probability distribution of 
all possible price paths between now and a certain future date. A useful model 
should therefore provide a well defined procedure to filter the series past price 
changes, and to compute the weights of the different future paths. Multifractal 
models offer a well defined set of answers to these questions, and their predictive 
power have been shown to be quite good |31[ 148 [ I32j . Other models, more in 
the spirit of GARCH or of Eq. (|15|) have been shown to also fare rather well 
|49| I41[ I44[ I47j . Filtering past information within both frameworks is actually 
quite similar; it would be interesting to compare their predictive power in more 
details, with special care for non-Gaussian effects. Once calibrated, these models 
can in principle be used for VaR estimates and option pricing. However, their 
mathematical complexity do not allow for explicit analytical solutions (except in 
some special cases, see e.g. [13]) and one has to resort to numerical, Monte-Carlo 
methods. The difficulty for stochastic volatility models or GARCH models in 
general is that the option price must be computed conditional to the past history 
|5Uj . which considerably complexifies Monte-Carlo methods, in particular for path 
dependent options, or when non-quadratic hedging objectives are considered. In 
other words, both the option price and the optimal hedge are no longer simple 
functions of the current price, but Junctionals of the whole price history. Finding 
operational ways to generalize, e.g. the method of Longstaff and Schwartz |51j . 
or the optimally hedged Monte-Carlo method of [^2j to account for this history 
dependence, seems to us a major challenge if one is to extract the best of these 
sophisticated volatility models. The toolkit that Mandelbrot gave us to describe 
power-law tails and long-range memory is still incomplete: is there an Ito lemma to 
deal elegantly with multifractal phenomena? Probably not - but as Paul Cootner 
pointed out long ago in his review of the famous cotton price paper Mandelbrot 
promised us with blood, sweat, toil and tears. After all, a faithful model, albeit 
difficult to work with, is certainly better than Black and Scholes's pristine, but 
often misleading, simplicity |53*j . 
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Figure 3: Multifractal scaling analysis of S&P 500 returns. 3-a: S&P 500 index 
price series sampled at a 5 minutes rate. 3-b: First five absolute moments of the 
index (defined in Eq. (j2j)) as a function of the time period r in double logarithmic 
representation. For each moment, a linear fit of the small scale behavior provides 
an estimate of C g . 3-c: Moment ratios M (J (t)/M2(t) <?//2 in log-log representation. 
Such curves would be flat for a geometric Brownian process. 3-d: Q spectrum 
estimate versus q. Negative q values are obtained using a wavelet method as 
defined in e.g. [2E]- This function is well fitted by a Kolmogorov log-normal 
spectrum (see text). 
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Figure 4: Multiplicative construction of a Mandelbrot cascade. One starts at the 
coarsest scale T and constructs the volatility fluctuations at fine scales using a 
recursive multiplication rule. The variables W are independent copies of the same 
random variable. 
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Figure 5: Multifractal properties of BMD model. 5-a: Exponential of a BMD 
process realization. The parameters have been adjusted in order to mimic the 
features of S&P 500 index reported in Fig. |21 5-b: Multifractal scaling of BMD 
return moments for q = 1,2,3,4,5. 5-c: Estimated ( q spectrum (o) compared 
with the log- normal analytical expression (solid line). 5-d: Evolution of the re- 
turn probability distributions across scales, from nearly Gaussian at coarse scale 
(bottom) to fat tailed law at small scales (top), (e) Log- volatility covariance as a 
function of the logarithm of the lag r. 
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vol_h 

Figure 6: The correlation between past historical volatilities cx/Jr] (ln(r) on the 
horizontal axis) and the future realized volatilities ovfr'] (ln(r') on the vertical 
axis) for the USD/CHF foreign exchange time series. The units correspond to the 
time intervals: lh 12min for (1), 4h 48min (4), Id 14h 24min (10), 6d 9h 36min 
(14), and 25d 14h 24min (18). 
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Figure 7: As for Fig. El but for a theoretical volatility cascade with Ornstein- 
Uhlenbeck partial log-volatilities £„. The cascade includes components from 0.25 
day to 4 days. Note the symmetry of this mug-shot, to be contrasted with the 
previous figure. A similar, symmetric mug-shot would also obtain for multifractal 
models.. 
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